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Ga interstitials in GaAs (Igo) are studied using the local-orbital ab-initio code SIESTA in a su- 
percell of 216+1 atoms. Starting from eight different initial configurations, we find five metastable 
structures: the two tetrahedral sites in addition to the 110-split[G a -As] , lll-split[Ga-As]j an d 
100-split[ Ga -Ga]- Studying the competition between various configuration and charges of Ic a at 
T = K, we find that predominant gallium interstitials in GaAs are charged +1, neutral or at 
most —1 depending on doping conditions and prefer to occupy the tetrahedral configuration where 
it is surrounded by Ga atoms. Our results are in excellent agreement with recent experimental 
results concerning the dominant charge of Igci, underlining the importance of finite size effects in 
the calculation of defects. 
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I. INTRODUCTION 

Gallium self-interstitials are believed to play a signifi- 
cant role for dopant diffusion in GaAs. The in-diffusion of 
an acceptor dopant A~\ (at an interstitial position) occurs 
via a kick-out mechanism that transforms it to a substi- 
tutional atom (AQ a ) and a gallium interstitial (lQ a ) plus 
the emission of a number of holes (equation from cited 
Ref.[Q): 
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where k denotes the charge state of J^ a involved in the 
reaction. 

Early calculations for Ga self-interstitials in GaAs^' 
led experimental groups to propose contradicting conclu- 
sions regarding the charge state of active self-interstitials 
in GaA s. The su ggested states var ied fr om neutraP, to 

+1 eeei + psno] or both +2 and + £nn] Recently? how _ 

ever, Bracht et a/P found that fits of recent as well as 
earlier diffusivity profiles are more accurate for dominant 
neutral and +1 charge states. This analysis of published 
data is confirmed by the observed compatibilities be- 
tween the hole concentration measurements and dopant 
(Zn) concentrations^. 

These experimental results demonstrate the need for 
a set of more detailed and accurate quantum mechan- 
ical calculations regarding the dominant charge state 
and geometry of Ic a m GaAs. Most recent papers only 
treat a sub group of the charge states and interstitial po- 
sitions^^, however, and we still lack a complete de- 
scription of the competition between different Ga self- 
interstitials in GaAs. This paper intends to fill this gap 
by providing a unified analysis of all charge states from 
q = — 3 to q = +3 for a wide range of the Ga self- 
interstitial configurations at T = K, in order to identify 
the dominant defects but also characterize others that 
could play a role in the diffusion of Ic a or after ion beam 
implantation, for example. 

This paper is organized as follows. Section \H\ explains 
the methodology used for defect calculation. Next, we 



describe in Section III the Gallium interstitial configura- 
tions used as starting points for this work. Section IV is 



devoted to study the stability of the chosen gallium in- 
terstitials after full relaxation of both the neutral and the 
charged states. The most relevant Ga interstitial config- 
urations and charge states in GaAs are then deduced and 
compared with earlier results in Section [V| 



II. METHODOLOGY 

A ll ca lculations are performed using the SIESTA 
cod d 15 * 16 -! within density functional theory (DFT) in local- 
density approximation (LDA). The details of the proce- 
dure followed is discussed in our previous worl^ and we 
focus below on the operations and parameters specific to 
the Ga self-interstitial simulations. 



A. Simulation parameters 

Simulations are performed using a supercell with 
216+1 atoms. This size is just sufficient to prevent size 
effects from dominating the structure and energetics of 
defects in GaAs. The wavefunctions are constructed from 
a double-^ polarized basis set (DZP) and we use a 2 x 2 x 2 
Monkhorst-Pack mesh sampling^. The choice of these 
parameters is discussed at length in our earlier work and 
the reader is referred to Ref . [T71 for more details. 

To test the validity of the local basis set used in this 
work, we evaluate the heat of formation of bulk GaAs 
crystal (Aif), defined as : 
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For this, it is necessary to compute the bulk chemical po- 
tentials, calculated from a 32 atoms As lattice a 
64 atoms Ga lattice and a 216 atoms GaAs lat- 

tice (fJaaAs) respectively. The heat of formation repre- 
sents the energy necessary to dissociate the GaAs crystal 
into its Ga and As components. Table [I] compares the 
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chemical potentials obtained using DZP with chemical 
potentials derived from a plane wave basis set calcula- 
tion (PW) within the DFT-LDA carried out by Zollo et 
a/P^on 64+1 atoms supercell. Our calculations provide 
an excellent agreement with experiment: both for the 
lattice parameter and the formation enthalpy. 



TABLE I: Comparison between DFT-LDA calculations — 
with double-^ polarized basis set (DZP) from this work and 
plane waves basis set (PW) from the work of Zollo et a/P^ 
— and experiment for the lattice parameter (a), chemical 
potentials (fi) and the resulting formation enthalpy (AH). 
Ab-initio calculations are performed at K and experimental 
data at 300 K. 





DZP 




ExptP 


a (A) 


5.60 


5.55 


5.65 


VGa (eV) 


-61.487 


-61.785 




VAs k (eV) 


-173.83 


-173.75 




^gTas (eV) 


-236.05 


-236.12 




AH (eV) 


-0.737 


-0.985 


-0.736 



B. Formation energy calculations 

Ga self-interstitials are placed at various sites in the 
supercell. Since these positions do not necessarily cor- 
respond to a local minimum, the network is slightly dis- 
torted and the configuration is relaxed at T = K until 
a total force threshold of 0.04 eV/A is reached. 

The formation energy (Ef) of each self-interstitial is 
calculated using 
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E f + q(E v + ii e ) - \(n As - n Ga )A/i (3) 



where Ef is the formation energy independent from dop- 
ing and growing conditions, the next term on the right- 
hand side depends on the doping of the sample \i e (i.e. 
Fermi level), the charge state of the defect q and the po- 
sition of the valence band maximum, Ey\ the last term 
is associated with the stoichiometry of the supercell con- 
taining riAs Arsenic and nca Gallium atoms. Finally, the 
chemical potential difference Afi is defined as: 
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The independent formation energy can thus be calcu- 
lated numerically using the relation: 



E f = E tot (q)--(n A s+n Ga )vGa A s 
1 
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where E tot (q) correspond to the total energy of the re- 
laxed supercell containing the self-interstitial. 



The total energy of the relaxed supercell must be cor- 
rected for the strong perturbation produced by the net 
charge on the relaxed state symmetry and local elec- 
tronic properties of the supercell. We can account for 
the electrostatic interaction between the charged defect 
and its periodic images by adding a neutralizing jellium 
background then correcting the relaxed energy (E tot (q)). 
Madelung correction due to the periodic boundary con- 
ditions is introduced following the Makov and Payne 
approximate procedure 20 . According to our previous 
WO rlp3 the monopole-monopole interactions correction 
is calculated to be 0.094 eV, 0.37 eV, and 0.84 eV for 
charge states ±1, ±2, and ±3 respectively, while higher 
order corrections were found to be negligible. Charged 
state formation energies, from Section |IV A| and later, 
were adjusted using these corrections. 

Finally, the position of the Fermi level fi e varies with 
doping and temperature and depends strongly on the 
carrier concentration. Thus, majority carriers (electrons 
or holes) can get trapped at defect levels changing the 
charge state of a given defect from qi to #2- The thermal 
ionization energy from a charge q\ to q^ is by definition 
the value of the Fermi level where the transition occurs: 
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We use Eq. [6] in Section IV C to calculate ionization en- 
ergies of charged defects for metastable configurations. 



III. INITIAL CONFIGURATIONS FOR THE GA 
INTERSTITIAL 

We first determine the metastable configurations for 
the Ga interstitial (Igo) with different charge states. We 
start in each of eight positions, relaxing the interstitial 
and characterizing the local energy minimum reached. 
All the initial states are shown on the top row of Fig. 
IT] The first starting point, from the left, tetra[Ga-As]? 
nas the Ga interstitial placed in a tetrahedral position 
with four surrounding lattice As atoms. For the second 
starting configuration, tetra[Q a -Ga] ? the lc a is shifted in a 
tetrahedral position with respect to 4 Ga atoms. In these 
two initial states, the bonds between the interstitial atom 
and its four tetrahedral neighbors, positioned on a perfect 
tetrahedron, have the same length of 2.425 A, which is 
exactly the length of Ga— As bonds in the zincblende 
structure corresponding to the lattice parameter we use 
a = 5.60 A (see Table 0. 

The third starting point is an hexagonal interstitial 
configuration (hexag), where Ia a is located at the cen- 
ter of the six-membered ring with alternating chemical 
species at equal distances from its six nearest neighbors. 
The Igcl— Ga and lea - As bond lengths are 2.322 A while 
^Ga — Ga— As angles all have the same value of 63.0°. We 
also examined a bond- center configuration in which lc a is 
lying exactly in the middle of a Ga— As bond at 1.213 A 
from each of them. 
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Tetra[Ga-As] Tetra[Ga-Ga] 



Hexag Bond-Center I I 0-SpNt[Ga-Ga] I 1 0-SpHt[Ga-As] I 00-SpNt[Ga-Ga] I 00-SpHt[Ga-As] 




E(Tetra[Ga-As]) = 3.73eV 



E(Tetra[Ga-Ga]) = 3.42eV E( I 1 0-split[Ga-As]) = 4.24eV E( 1 00-split[Ga-Ga]) = 4.92eV 



FIG. 1: (Color online) Top: The eight initial configurations considered in this study for neutral self-interstitial Igu- The six 
first configurations, going from left to right, are viewed near the <110> direction, while the remaining two are viewed along the 
<100> direction. Bottom: The metastable configurations obtained after full relaxation of Ic a - Full arrows connect the initial 
configuration to its metastable counterpart, while dashed arrows means that the initial configuration is unstable and converged 
to the pointed configuration. Gallium atoms are red while arsenics are white; the interstitial Ga atom is colored yellow. For 
splits, regular sites of the displaced lattice atoms are highlighted by a dotted circle. 



Finally, we look at four different interstitials from the 
important family of split geometries. Split interstitials 
are formed when Ic a pushes one regular lattice atom (Ga 
or As) out of its crystalline position, forming a dumbbell 
centered at a regular lattice site. The split-interstitial 
type is determined by the orientation of the vector join- 
ing the pair of atoms (see Ref. 21). Gallium interstitials 
can form dumbbells with As and Ga atoms following the 
<100> and <110> directions, as shown in the top part 
of Fig. [I] or the <111> direction. 

In more details, in the 110-split[G a -As] interstitial, an 
arsenic atom is moved by 1.570 A from its regular lattice 
site to make room for the intserstitial atom positioned at 
0.840 A from the regular lattice site and forming a dumb- 
bell length of about 2.312 A along the <110> direction. 
In this case, the center the dumbbell is slightly displaced 
from regular lattice site (as clearly seen in Fig.|2ja)). 

For 110-split[Ga-Ga] 7 a lattice Ga atom is moved by 
1.338 A from its lattice site along the dumbbell axis (in 
another <110> direction) and the interstitial is placed 
at 0.641 A in the opposite direction along this axis from 
the lattice site. The dumbbell length is now 1.980 A. 



For 100-split[Ga-Ga] an d 100-split[Ga-As] interstitials, the 
lattice atom (Ga and As respectively) is moved along 
the dumbbell axis from its regular position by 1.212 A 
while the interstitial is placed at the same distance in 
the opposite direction along the same axis, forming a 
dumbbell of 2.425 A. 

Other possible configurations including interstitial 
clusters might exist in real crystals. As a first step, we 
restricted ourselves to these simplest structures. 



IV. RESULTS 

Here we present the results of our simulations, us- 
ing the techniques and parameters described in Sec- 
tion [IlJ We first discuss the stability of the eight inter- 
stitial positions described in Section |III| in the neutral 
state. Then, we focus on the influence of the charge 
state on the metastable interstitial positions identified 
earlier. Finally, we discuss the impact of the progressive 
doping of the material on the competition between vari- 
ous charge states of a given interstitial in stoichiometric 
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GaAs (A/i = ). 



A. Structural stability of neutral self- int erst it ials 

The bottom part of Fig. [I] shows the final geometry of 
the relaxed interstitial configurations in the neutral state. 
The most stable interstitial is the tetra[G a -Ga]? which 
undergoes small lattice distortions leading to its conver- 
gence into a metastable configuration with a formation 
energy of approximately 3.42 eV. The hexag, bond-center 
and 110-split[G a -Ga] configurations are unstable and re- 
lax to the same tetra[Q a -Ga] after undergoing consider- 
able atomic displacement and lattice relaxation. 

The second tetrahedral configuration, tetra[Ga-As]? 
is also metastable, with a formation energy of about 
3.73 eV, slightly above that of tetra[G a -As]- This struc- 
ture is only a few relaxation steps away from the initial 
tetra[G a -As]- Both tetrahedral interstitials leave the sur- 
rounding crystalline network relatively unaffected. The 
tetra[Ga-Ga] configuration is close to the starting config- 
uration with only a slight outward relaxation of the sur- 
rounding Ga neighbours leading to an increase in length 
of both Ica-Ga bond by about 7.0 %, to 2.596 A. The 
volume expansion around the tetra[G a -Ga] g° es down 
rapidly and affects only the first and second neighbor 
shells along the tetrahedral axes. The tetra[G a -As] con_ 
figuration undergoes a similar expansion and the Igo - As 
bond lengths are stretched by 5.3 % from 2.425 to 
2.554 A. 

Both the 110-Split[G a -As] an d 100-Split[G a -Ga] are 
found to be metastable. The formation energies are 
higher: E f = 4.24 eV for 110-split[ Ga -As] and E 'f = 
4.92 eV for 100-split[G a -Ga]- I n addition, the stress im- 
posed on the lattice is more important, and affects sig- 
nificantly the more distant neighbours. 

The 100-split[G a -Ga] experiences the largest lattice de- 
formation around the defect among other interstitial de- 
fects as illustrated in Figure [I] where surrounding lattice 
atoms are shown. This relaxed configuration is a dumb- 
bell formed by Ic a an d the displaced Ga lattice atom. 
The dumbbell is centered and symmetric with respect to 
the middle vacant lattice gallium site with each Ga atom 
being located at 1.10 A apart. The length of the dumb- 
bell shrinks by about 9.2% (from 2.42 to 2.20 A) bringing 
the two atoms closer. For their part, atoms at the first 
and second shell neighbors experience an outward relax- 
ation and are pushed away from their original position 
by about 0.46 and 0.21 A respectively. Thus, consid- 
ering both effects the bond length between each of the 
Ga atoms forming the dumbbell and their first As lattice 
neighbors increases by approximately 17.4% from 1.99 to 
2.34 A as the distance to the second nearest neighbors 
grows from 3.22 to 3.47 A, a change of about 7.8%. 

Finally, the 100-split[G<a-As] self-interstitial is highly 
unstable and does not converge to any stable state. For 
this reason, we will not attempt any further calculation 
using this configuration for the rest of this work. 



TABLE II: Formation energies (in eV) for stable and 
metastable Ga interstitial configurations in GaAs for various 
charge states. 



Stable 






Net system charge 


Q 




configurations 


-3 


-2 


-1 


+ 1 


+2 


+3 


tetra [G a-Ga] 


8.37 


6.40 


4.62 3.42 


2.40 


2.81 


3.43 


tetra [Ga _ As] 


8.65 


6.72 


4.98 3.73 


2.67 


2.91 


3.41 


110-split [Ga _ As] 




6.42 


4.99 4.24 


1 


1 


1 


lll-split [Ga _ As] 








3.33 


3.74 


4.35 


100-Split[Ga-Ga] 




7.75 


6.16 4.92 


3.33 


4.15 


4.62 



In order to test the stability of the four metastable con- 
figurations found, we further relaxed them with a more 
accurate force threshold of 0.002 eV/A. The observed 
change in geometry and formation energy is negligible in- 
dicating that our results are already well converged with 
the former force threshold of 0.04 eV/A . 



B. Structural stability of charged self- interstitials 

Having characterized the stability of the set of neu- 
tral initial self-interstitial configuration, we now turn to 
charge states. All charge configurations are also started 
from the eight ideal unrelaxed configurations except for 
the 110-split[G a -As]? which could not relax in the allowed 
time from the ideal position and which was started in the 
neutral relaxed configuration instead. 

For q = +1 charged interstitials, the relaxation fol- 
lows the same scenario as for the neutral defect for all 
first seven configurations of Fig. [T] (top): the overall 
stability order is kept unchanged and the tetra^ a _ Ga j 
configuration is still the most stable defect. For their 
part, the unstable interstitial states relaxed into the same 
metastable configuration as in the neutral case with one 
exception, however: the 110-split[G a -As] now relaxes into 
a Hl-split[G a -As] (Fig. [5| but keeps the same stability 
order with respect to the other metastable defects. 

Interestingly, the resulting formation energies calcu- 
lated for +1 charged defects are significantly lower than 
in the neutral charge state in all cases: tetra^ a _ Ga j 

has a formation energy of Ej = 2.40 eV, followed by 

tetra+^^j with E f = 2.63 eV, then lll-split+ 1 a _ As] 

with Ef = 3.33 eV and finally lOO-split^^ with also 

E' f = 3.33 eV. 

The removal of an electron for lQ a in GaAs stabilizes 
all defects uniformly with respect to their respective neu- 
tral state while the stability order of each interstitial 
configuration with regard to each other remains about 
the same. To fully characterize this effect, we have fur- 
ther relaxed the most stable interstitial geometries — 
tetra [Ga -Ga], tetra [Ga -As], 100-split [Ga -Ga], and 110 or 
111 splitS[G a -As] — f° r the q = — 1, q = ±2 and ±3 
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FIG. 2: (Color online) The two distinct 110 (a) and 111 (b) 
splits seen from a 110 view. Dash lines correspond to <111> 
and <110> crystalline axis (indicated in blue) and dotted 
circles refer to the regular crystalline position of the displaced 
arsenic atom. Gallium atoms are red while arsenics are white; 
the interstitial Ga atom is colored yellow. Some lattice atoms 
have been removed for clarity. 



charges, supposing that the unstable interstitial configu- 
rations do not stabilize in these highly charged systems. 

Table |TT] shows the formation energies of the five 
metastable relaxed configurations in increasing order 
of formation energy, for seven charge states (q = 
±3, ±2, ±, and 0) with the associated monopole correc- 
tion applied, as explained in Section |IIB| The missing 
numbers for charge q = — 3 correspond to configura- 
tions that did not achieve convergence even after long 
simulations. Additionally, starting from the the neutral 
110-split[G a -As] configuration, relaxations towards pos- 
itive charges all induce a change in the orientation of 
the dumbbell leading to the nearly same Hl-split[G a -As] 
(as indicated by downside arrows in Table |n| while re- 
laxations of negatively charged 110-split[G a -As] preserve 
the initial <110> orientation. 

Figure [2] illustrates the shift of orientation from 

110- split[Ga-As] f° r neutral and negative charges to 

111- split[G a -As] f° r positive charges. The length of the 
dumbbell is around 2.24 A for lll-split [G a-As] +1 ' +2 ' +3 
and about 2.33 A for 110-splitfGa-As] '" 1 '" 2 - 

From Table [TTJ we observe that the lowest formation 
energy for all the four interstitial configurations is asso- 
ciated with the +1 charged state. For a given symmetry, 
we see that the formation energy monotonically decreases 
with increasing charge, from —3 to +1 before going up 
for more positive charge states. In all cases, the forma- 
tion energies for positive charges are lower than for the 
negative ones. 

Looking at particular charged interstitials, we ob- 



serve also that the order of stability between differ- 
ent interstitial configurations that we observe in neu- 
tral charge state is conserved under the variation of the 
net charge of the system (except for —2 charged state 
where 110-split[G a -As] formation energy is below that of 
tetra[Ga-As])- Moreover, apart from higher charge states, 
there is a somewhat constant difference of 0.30 eV be- 
tween the two different tetrahedral configurations for the 
same charge states from —3 to +1 where the formation 
energy reaches its minimum value. 

We note also that a degeneracy appears between pairs 
of interstitials for two charge states: at +1, both 111 and 
100 splits have about the same formation energy while at 
— 1, the same is true for tetra[G a -As] an d 110-split[G a -As]- 
This suggests that the 110 and 111 split [Q a - As] configu- 
rations play a key role in GaAs crystals, serving as tran- 
sitional configurations when passing from tetra[G a -As] a ^ 
charge —1 to 100-split[G a -Ga] at charge +1, successively 
losing two electrons, one at a time (the Ill-split being 
an intermediate step from neutral to +1 charged state). 
As a result, this specific transition process could be an 
important diffusion path for impurities in GaAs crystals. 

Although, for charge +3, the formation energies of 
both tetrahedral configurations are also near, it could 
be a finite size effect associated with the stress induced 
by the high positive net charge, introducing a bias in our 
calculation. 



C. Competition between Ic a charge states under 
doping conditions 

We now look at the effect of doping by varying the 
Fermi level with the help of the parameter ja e in equa- 
tion [3J These effects are best seen by comparing data for 
multiple charge states and we concentrate on the con- 
figurations of Table [TTJ Because of the similarities in the 
stability diagrams of tetra[G a -Ga] an d tetra[G a -As] on ly 
the first one is shown in Figure [3] The dia gram for 
100-split[Ga-Ga] an d 1 10-split [G a - As] are shown in Fig- 
ures [4] and [5] respectively. In all figures, the Fermi level 
is set by reference to the valence band maximum. 

Because DFT calculations are known to underestimate 
the band gap (the present calculation gives a band gap 
of 0.82 eV compared to 1.52 eV at T = K reported 
by experiment), it is common to vary the Fermi level in 
the window of the experimental band gap to obtain the 
full picture^. Note, that even if our calculated ioniza- 
tion levels have a conduction band character, they have 
not been corrected for that because the efficiency of the 
several band gap correction methods is system depen- 
dent and testing them is out of the scope of the present 
wor ]p9 

Fig. [3] shows the formation energy as a function of 
the doping level for the tetrahedral insterstitial. Dot- 
ted lines point at the location of ionization level, iden- 
tified by intersecting formation energy lines. Three sta- 
bility domains, labeled by (I), (II), and (III), are found 
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FIG. 3: (Color online) Formation energies as function of Fermi 
level for various charge states of the tetra[G a -Ga] configuration 
at K. Dotted lines point at the ionization level locations 
delimiting the three stability domains labeled by (I), (II), and 
(III). 



FIG. 4: (Color online) Formation energies as function of Fermi 
levels for various charge states of the 100-split[ Ga -Ga] config- 
uration at K. 



for tetra[G a -Ga] corresponding to successive dominant 
charge states +1, 0, and —1 respectively. The same sta- 
ble states occur for tetra[£ a _^ s ] , with almost the same 
ionization energies. Table |III| summarizes these ioniza- 
tion energies, calculated from equation [6j 

The 100-split[G a -Ga] exhibits a different behavior as 
the line of charge state q = — 1 crosses the line from the 
state q = +1 before the horizontal neutral q = state 
line (about 0.18 eV below) allowing a direct transition 
between q = — 1 and +1, in a so-called negative- U effect. 
Only two stability domains labeled by (I) and (II) are 
found, for charge +1 and —1 respectively (see Fig. [4]). 
This effect might not exist in real systems since region II 
occurs at the edge of the conduction band. Any fluctua- 
tion or error might therefore screen region II. The reverse 
could be true for the tetrahedral interstitials (Fig. [3]). In 
this figure, region II is very narrow, set in the middle 
of the gap, and a fluctuation could remove it altogether, 
this time inducing a negative-U effect. 

We note also that the transition domain of the 100- 
split occurs just below the minimum of the conduction 
band, meaning that it is only accessible in extreme doping 
condition. For most purposes, the domain of charge +1 
will be the only one that matters. 

Another negative-U effect causes the transition from 
lll-split [Ga _As] +1 to 110-spHtjp a _ As] (Fig. [5]). Once 
again, the transition occurs very near to tne neutral 
charge state line, only 0.07 eV below. Here, however, 
the transition is located at midgap and should therefore 
play a more important role. Moreover, as this (+1/— 1) 
negative-U transition manifests itself by a change in the 



TABLE III: Ionization energies of metastable Ioa configu- 
rations in GaAs (see equation [6]). Negative-U transition 
(+1/ — 1) for split [ Ga _ As] changes the dumbbell orientation 
from 111 to 110. 







Ionization levels (eV) 




Configurations 


+1/0 


0/-1 


Negative- U 
+ V-1 - 


-1/-2 


tetra [Ga _ G a] 
tetra [Ga -As] 


1.02 
1.06 


1.20 
1.25 






100-Split [Ga -Ga] 






1.42 




110 & 111 splits [Ga _As] 






0.83 


1.42 



orientation of the dumbbell from 111 to 110, it should 
therefore be relatively insensitive to the various limita- 
tions of the current simulation and other possible thermal 
fluctuations. 

Similarly to the 100-split, we find also a transition 
(here, — 1/— 2) for 110-split[G a -As] near the conduction 
band minimum. Since this transition level has itself 
a conduction band character its position might be af- 
fected by the DFT-bandgap underestimation and must 
be treated with care. Although we did not manage to 
converge the 110-split[G a -As] -3 ? it is interesting to note 
that if the trend for the negative charges seen in Fig. [5] 
holds, we should see a negative-U transition between 
charge —1 and —3 before the — 1/— 2 transition for this 
charge. 
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FIG. 5: (Color online) Formation energies as function of Fermi 
levels for various charge states of the 110 (q = —2,-1,0) and 
111 (q = +l,+2,+3) split[Ga-As] configuration at K. 



V. DISCUSSION 

Here we discuss and compare our res ults with p revious 
ab initio and tight-binding calculation d 13 * 14 * 1 -^^, as well 
as with recent experimental data from BrachtP'. While 
it was not clearly indicated in most of the theoretical 
works, whether the tetrahedral interstitial label meant 
tetra[G a -As] or tetra[Q a -Ga] ? we will presume that it refers 
to the former. 

Chadi 21 was first to report self-interstitial config- 
urations and energetics using self-consistent pseudo- 
potential relaxations on GaAs supercells with 32+1 
atoms. Almost the same set of starting configurations 
as in the present work was studied under different charg- 
ing, but the resulting stability order was completely dif- 
ferent from our. Indeed, Chadi found the bond-center 
(twofold) configurations to be the most stable configura- 
tion for q = +1 and the 110-split[G a -As] f° r Q = 0,-1. 
The tetra[Qa-As] was found to have the lowest formation 
energy only under +2 charging. The difference between 
this work and ours is mostly due to the strong finite size 
effects present in a 32+1 unit cell. 

For their part, Zollo and NiemineiP^ studied the full 
set of interstitial positions — except for tetra^Qa-Ga] — m 
a 64+1 atomic supercell and for the neutral state only. 
They found that the hexagonal, bond-center, and 100- 
splits are unstable, converging to the tetrahedral inter- 
stitial position. They also identified the tetra[Ga-As]? 
the 110-split[G a -Ga]? an d the 110-split[G a -As] to be 
metastable, with increasing formation energy. Recently, 
Volpe et a/P^l used a large supercell of 216+1 atoms 
with tight binding method to study lea i n GaAs, again 
treating exclusively the neutral charge state and com- 



puting formation energies relative to the tetrahedral in- 
terstitial formation energy only. The metastable struc- 
tures identified at the neutral state were classified in 
increasing order of formation energy: the tetra[G a -As] ? 
then the 110-split[G a -Ga] •> the 110-split[G a -As] > an d the 
100-split[G a -Ga]- As was shown in Tab. |TTJ our results 
show that 100-split[G a -Ga] is indeed metastable as found 
by Volpe et al. They disagree with both Zollo and Niem- 
inen and Volpe et al. with respect to the 110-split[G a -Ga] 
interstitial, however, which becomes unstable and prefers 
to relax to tetra[G a -As] according to our 216 DFT calcu- 
lations. 

The difference between these calculations and the one 
presented here are caused by (1) size effect associated self- 
interactions of the defects in unit cells that are too small 
and (2) the accuracy of the potential (DFT versus tight- 
binding), particularly when important structural changes 
are taking place, as is the case for split interstitials, for 
example. 

Focusing on the dominant charge state of Ia a in GaAs, 



the results of Section [TV] show that higher charge states 
(q = ±2, ±3) are not relevant and should contribute neg- 
ligibly to the total experimental self-diffusion profiles. 
This agrees well with recent experimental results from 
Bracht 1 that identified J^ a and IqI as important species 
for diffusion process in GaAs crystals doped with Zn (but 
with lower contribution than vacancies) . Our results pro- 
vide a strong support for the picture proposed by Bracht, 
disproving earlier models which generally predict diffu- 
sion process controlled mainly by +2 and +3 intersti- 
tials. In particular, our calculations show that +2 and 
+3 charge states exhibit higher formation energies than 
the +1 charge defect, contrary to what was found by 
Zhang and NorthurrM Using a 32+1 atoms supercell 
within DFT-LDA, these authors identified the dominant 
native defect to be the tetrad 3 A , in Ga-rich condition 

[Ga- AsJ 

under p-type doping. Similarly, more accurate ab-initio 
calculation of Zollo et al. 18 find negative-U effects for 
tetra[G a -As] with (+3/+1) and (+1/— 1) ionization levels 
located at 0.29 and 1.23 eV above the VBM respectively, 
still giving an important role to the triply positive state 
which we do not see. 

Again, size effects can explain many of these diver- 
gences. For example, finite size effects have been reported 
recently by Schick et al. 24 for As 110-split interstitials in 
GaAs. While a 65 supercell calculation with a 2 x 2 x 2 
k-points mesh predict the stability of +2 charge state 
starting for the VBM. This charge state disappears com- 
pletely from the diagram as soon as a supercell as large 
as 217 is used with different k-point meshes leaving the 
+1 charge state as the most stable near VBM. 

Qualitatively, the formation energies for tetra[G a -Ga] 
in stoichiometric GaAs (A/i = 0) we compute depend on 
doping conditions: E f (I% a ) = 3.42 eV, E f (I+l) = 2.4 
to 3.92 eV, E f (I+ 2 a ) = 2.81 to 4.33 eV, and for 

E f( I Ga) = 3 - 43 to 4 - 95 eV - A11 tnese values remain in 
the window of allowed values compared to the activa- 
tion enthalpy (H a ) obtained after fitting the experimen- 
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tal profiles. H a for / G 'J~ -mediated Ga diffusion in GaAs 
reported by Bracht 1 was 5.45±0.12 and 5.80±0.32 eV for 
neutral and +1 charge states respectively. Since we do 
not know the migration enthalpies of Io a in GaAs with 
respect to the charge state, it is not possible at this point 
to push further and identify the charge state responsible 
for the Zn diffusion profiles. Only a detailed study of the 
migration mechanisms of Ia a in GaAs similar to the one 
performed recently for Vca m GaAs as function of the 
charge state 25 can give the answer to this question. 

VI. CONCLUSIONS 

In this work, we have studied the stability of 
Ga self-interstitial for multiple charge states (q = 
0,±l,±2,and ± 3) within DFT-LDA using the local- 
orbital basis set program SIESTA at T = K. Out of 
the eight initial configurations tested, five were found to 
be metastable after full relaxation. As a general rule, 
the most stable configuration is found to be tetra[Q a -Ga] 
for all charge states, in addition positively charged inter- 
stitials are more stable than negative ones for all tested 
charges and configurations. 

After studying the competition between various con- 
figuration and charges of lea, we conclude that predom- 
inant gallium interstitials in GaAs are charged +1, neu- 
tral or at most —1 depending on doping conditions. This 
agrees well with the recent conclusions driven by Bracht 
et at. which states that fits of recent as well as ear- 
lier diffusivity profiles in Zn doped GaAs are more ac- 
curate if the role Iq^° is considered. At low tempera- 
tures, when the formation energy dominates, lQ a prefers 
to occupy the tetrahedral interstitial configuration being 
surrounded by gallium atoms (tetra[Q a -Ga])- The com- 
petition between tetra[Q a -Ga] an d the other metastable 
configuration increases as we approach the experimental 



processing temperatures (above 1000 K) but not suffi- 
ciently to invert the order of stability. For example un- 
der P-type doping at 1000 K, the tetra[Q a -Ga] +1 still has 
a probability of occurrence about 100 times larger than 
and 10 5 times larger than the 110 and 100-split in- 
terstitials. 

The comparison of our results with previous works 
show also that size of the simulation supercell can affect 
significantly the stability and formation energy of lea, 
and both tight-binding and ab-initio calculations become 
more reliable with increasing cell size. The size of the su- 
percell affects the charge state of the dominant defect 
and modifies also the metastability of other defects such 
as the split interstitials. In this case, the use of a 216+1 
atoms supercell allows us to observe for the first time 
a change in the orientation of the split from <110> to 
<111> after the removal of one electron from the neutral 
configuration. 

In spite of the excellent agreement with recent exper- 
imental results, further calculations, including the en- 
tropic contributions and the diffusion pathways, are still 
necessary to obtain a complete picture of the role of Ic a 
and reveal the possible importance of other charge states 
and configurations. 
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